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Abstract 


The choice of high-lift system is crucial in the preliminary design process of a 
subsonic civil transport aircraft. Its purpose is to increase the allowable aircraft weight or 
decrease the aircraft's wing area for a given takeoff and landing performance. However, 
the implementation of a high-lift system into a design must be done carefully, for it can 
improve the aerodynamic performance of an aircraft but may also drastically increase the 
aircraft empty weight. If designed properly, a high-lift system can improve the cost 
effectiveness of an aircraft by increasing the payload weight for a given takeoff and 
landing performance. This is why the design methodology for a high-lift system should 
incorporate aerodynamic performance, weight, and cost. 

The airframe industry has experienced rapid technological growth in recent years 
which has led to significant advances in high-lift systems. For this reason many existing 
design methodologies have become obsolete since they are based on outdated low 
Reynolds number wind-tunnel data and can no longer accurately predict the aerodynamic 
characteristics or weight of current multi-element wings. Therefore, a new design 
methodology has been created that reflects current aerodynamic, weight, and cost data 
and provides enough flexibility to allow incorporation of new data when it becomes 
available. 
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1.0 Introduction 

The design of an efficient high-lift system remains as challenging today as it was 
twenty years ago when A.M.O. Smith wrote his enlightening papers on high-lift 
aerodynamics. 1 - 2 Modern civil transport aircraft require complex multi-element high-lift 
systems to meet stringent performance criteria during the takeoff and landing phases of 
flight. In the current competitive market place, new aircraft designs are driven to simpler, 
more efficient high-lift systems that provide improved aerodynamic performance in terms 
of increased maximum lift coefficient, C, , increased lift-to-drag ratio, L/D, or 

max 

increased lift coefficient, C L , for a given angle of attack and flap setting. Gamer and his 
co-workers at Boeing 3 present excellent examples illustrating the importance of a high-lift 
system in the design of a B777 type of aircraft: 

1 . AC l =+ 0. 1 0 for a constant angle of attack on approach for landing reduces the pitch 
attitude angle by about one degree. For a given aircraft geometry and landing-gear 
location this allows a reduction in landing-gear height and an associated weight 
reduction of 1,400 lb. 

2. AC l = + 1.5 % at a fixed approach speed results in an increase in payload of 
6,600 lb. 

3. A(L/D)=+0. 1 0% on takeoff results in an increase in payload of 2,800 lb. 

These examples illustrate the enormous importance of a well designed and 
engineered high-lift system in the overall development process of a subsonic civil 
transport aircraft. However, these systems also increase the structural weight, 
complexity, maintenance requirements, and cost of an aircraft. Thus, the designer is 
faced with the task of developing such a high-lifit system that allows the airplane to meet 
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the takeoff and landing performance requirements while minimizing the weight and cost 
of the airplane. 

This report presents a design methodology for multi-element high-lift systems for 
subsonic civil transport aircraft that includes aerodynamic performance, structural weight 
considerations, system complexity, and cost. The methodology is designed to be 
compatible with ACSYNT 4 , a multidisciplinary computer-aided aircraft conceptual 
design tool. Consequently, important considerations are (1) to find a reasonable 
compromise between the CPU requirements (on current generation workstations) of the 
high-lift module and the accuracy of the predictions by the module and (2) to provide 
enough flexibility to enhance the capabilities of the high-lift module when more powerful 
hardware and software becomes available. 

It is envisioned that this methodology will be integrated into ACSYNT as 
illustrated in figure 1.1. Once the development of the aerodynamic database for double 
and triple-slotted flaps is finished, the high-lift module can be used in the I s ' level of 
design as described in this report. Given the initial concept and mission requirements, 
iteration is required to find the optimum initial high-lift configuration based on general, 
or historical data. Subsequent levtls of design should then take advantage of 
computational fluid dynamics (CFD), experimental fluid dynamics (EFD), or flight data 
to more accurately model the configuration selected from the 1 st level. When provided 
with more complete aircraft geometry and given data that is based on a specific 
configuration, the methodology can then converge on an optimum preliminary high-lift 


configuration. 
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Figure 1.1 Overview of High-Lift Module Integration 
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2.0 High-Lift Configurations 

High-lift systems consist of leading and trailing-edge devices. Leading-edge 
devices increase the maximum lift of an airfoil by delaying its stall angle, Aa max , as 
shown in figure 2. 1 . This change in the stall angle is relatively constant for airfoils both 
with and without trailing-edge devices. 5 Trailing-edge devices produce a lift increment, 
AC l , and as illustrated in figure 2.1 the magnitude of this increment is approximately 
independent of a leading-edge device. 

2.1 Trailing-Edge Devices 

A trailing-edge device generates additional lift through an increase in the effective 

rearward camber of an airfoil. It is desirable to create a flap that produces a large increase 

in lift while maintaining a high lift-to-drag ratio, L/D, in order to enhance both the takeoff 

and landing performance of an aircraft. There are many trailing-edge devices employed 

today, but the most common are the plain flap, slotted flap, and Fowler flap. Figure 2.2 

* 

shows examples of these various flap configurations. 

In addition to the various flap configurations, there are three widely used support 
types: hinge, linkage, and track supports. These systems, which are illustrated in figure 
2.3, are used on the various configurations to deploy the flaps to their proper deflection 
angles and provide necessary structural support. The hinge support is the simplest of 
these support types, but it also has the worst aerodynamic performance. Since the flap 
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kinematics are restricted by the hinge position, the flap can only be optimized at one 
deflection angle. The linkage and track supports are more complex, but they can provide 
optimum flap settings at multiple deflection angles. 

2.1.1 Plain Flaps 

The plain flap consists of a hinged trailing-edge with a gap that is usually sealed 
to reduce leakage. It is the simplest of all trailing-edge devices but is also the least 
efficient in terms of the lift increment it can generate without flow separation. Plain flaps 
are commonly used as high-lift devices on general aviation aircraft which have low wing 
loadings. On civil transport airplanes plain flaps are rarely used except as control 
surfaces. 

2.1.2 Slotted Flaps 

Slotted flaps offer improved efficiency over plain flaps due to a delay in flow 
separation. To better understand this phenomenon, we must review A.M.O. Smith’s 
classic papers on high-lift aerodynamics 1-2 . Previously, the effects of slots in multi- 
element airfoils were incorrectly attributed to boundary-layer control through blowing. 
Smith argues that the principal effect of a slot is to delay flow separation through inviscid 
interactions. 

Smith asserts that gaps have five primary effects on multi-element airfoils. The 
peak pressure on a downstream element is reduced by an induced velocity created by the 
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circulation of a forward element. This is referred to as the “ slat effect ”, and it delays 
flow separation by relieving pressure recovery. In turn, the “ circulation effect” induces a 
greater circulation on the forward element because its trailing-edge lies in a region of 
high velocity at the leading-edge of the adjacent downstream element. Since the trailing- 
edge of the forward element lies in a high velocity region, the boundary layer is shed at a 
high velocity. This “ dumping effect ” delays flow separation by reducing the pressure 
rise over the airfoil and allows for “ off-the-surface pressure recovery ” which is much 
more efficient than recovery in contact with a wall. Finally, each surface benefits from a 
'fresh boundary-layer ” which originates from the leading-edge of each element. Thus, 
multi-element airfoils allow for a more efficient pressure recovery, since thin boundary 
layers can withstand stronger adverse pressure gradients than thick ones before 
separating. 

While these five effects primarily influence the inviscid nature of the flow, they 
also have important secondary viscous effects. It is this balance between the inviscid and 
viscous nature of multi-element airfoils that necessitates the optimization of gap size. 

The inviscid effects favor smaller gaps, while the viscous effects require larger slots. One 
problem associated with this trade-off is found in confluent boundary layers. If the gaps 
are not designed properly the wakes will merge resulting in exceptionally thick boundary 


layers. 
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2.1.3 Fowler Flaps 


Improved performance can also be obtained if a flap creates Fowler action, or 
rearward translation of the flap. Here, Fowler action is defined as the measure of the 
change in position of the leading-edge of the flap in the plane of the chord of the fore 
element. This is illustrated in figure 2.4 and is expressed as: 

Fowler action = Ac = s, + s 2 (1) 

The extended chord of an airfoil can then be defined as the length of the cruise airfoil plus 
the Fowler action: 


c' = c + Ac (2) 

This extension in wing area increases the airfoil's lift curve slope, generating more lift 
without a significant increase in drag: 


C, = C, - 

a a clean c 


( 3 ) 


However, Fowler action also produces an increase in the nose-down pitching moment 
which makes the aircraft more difficult to trim. 


Of all slotted trailing-edge devices, the single-slotted Fowler flap is the simplest 
and most efficient. It is a superior device for takeoff because it has the best L/D, but 
designers often cannot employ wings with single-slotted flaps because their usable lift is 
inadequate for takeoff and/or landing. In such cases a multi-element flap is required. 

The usable lift problem faced in the design of civil transport airplanes is depicted 
for the landing case in figure 2.5. The landing distance of an airplane is governed by its 
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weight and approach speed, assuming that the wing area is fixed and determined by 
cruise requirements. The Federal Airworthiness Requirements stipulate that the 
minimum approach speed V A equals 1 .22V S ^, where V s ^ represents the 1-g stall speed of 
the airplane in the landing configuration. For a given weight and landing distance the 
approach speed is more-or-less fixed, and the lift coefficient during landing approach is: 

c La = w L /(2pV A ! s) < c Lm J(m) 2 (4) 

In figure 2.5a this lift coefficient is marked. During approach, airplane angle of attack is 
governed by the following equation for steady flight: 

a = 9-y (5) 

where a typical glide slope angle is y =-3°. For many configurations the maximum 
attitude angle is severely limited by the tail scrape angle (Fig. 2.5b) and, consequently 
“limit = ^taii-scrape + 3° (Fig. 2.5a). Hence, the combination of a given approach speed and a 
limited angle of attack may force the designer to select a more complex high-lift system 
as shown in figure 2.5a. 

Further delay of flow separation at higher flap deflections, an additional increase 
in effective camber, and a potential for greater Fowler action are produced with additional 
flap elements. With all of these benefits combined a multi-element flap is able to 
generate a higher maximum lift. It also has a larger usable lift at lower angles of attack 
which may be important if the fuselage angle is restricted on takeoff or landing. But this 
increase in lift comes with a price, for the efficiency, or L/D, is reduced with each 
additional flap element. Structural complexity also increases with the number of 
elements, so multi-element flaps become more costly to manufacture and maintain 
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2.2 Leading-Edge Devices 

Leading-edge devices are used primarily to delay the onset of stall by reducing the 
peak velocity in the leading-edge region of the main element. This corresponds to an 
increase in C, due to a shift in the stall angle. Slats and Kruger flaps are the most 

‘max ^ 

widely used leading-edge devices in industry today. The various available configurations 
are illustrated in figure 2.6. 

There is some disagreement on which configuration provides the highest 
maximum lift or the lowest drag. Wedderspoon 6 cites that a vented slat was chosen over 
a sealed folding bullnose Kruger for the Airbus 320 because it generated a higher 
maximum lift, but he does not mention if a vented Kruger flap was investigated. 
Woodward and Lean 5 argue that vented Kruger flaps produce considerably higher 
maximum lift than vented slats and that sealed Krugers produce higher maximum lift than 
sealed slats. Such inconsistencies in opinions indicate that these devices are not fully 
understood. Clearly, more research must still be conducted to better understand the 
aerodynamics of leading-edge devices. It appears that slats and Kruger flaps have similar 
aerodynamic performance if properly designed. 

In this paper, Kruger flaps and slats are considered to have the same lifting 
effectiveness but slightly different maximum lift capabilities. For the purposes of this 
methodology, the focus is placed on the flap position rather than on the configuration 
type. Slotted, or vented, leading-edge devices have high maximum lift capabilities. As a 
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result these devices produce higher drag, so they are best utilized during landing. 
Conversely, sealed leading-edge devices have lower maximum lift which produce less 
drag and should be used during the takeoff phase. 

2.2.1 Slats 


A slat offers more flexibility than a Kruger flap. It generally has three settings: 
stowed, takeoff, and landing. The tracks are constructed to optimize the configuration for 
each maneuver, so the slat is sealed for takeoff and vented for landing. 

2.2.2 Kruger Flaps 

There are several types of Kruger flaps available: simple Krugers, folding 
bullnose Krugers, and variable camber Krugers, and each configuration can be either 
sealed or vented. The primary difference between this leading-edge device and a slat is 
the support structure. Kruger flaps are hinged and stow on the lower surface of the 
airfoil, while slats deploy from the' leading-edge of the airfoil on tracks. The Kruger’s 
hinged support does not allow a multitude of position settings, so it must be either 
retracted or extended. It is up to the designer to optimize the flap for takeoff or landing. 
Generally, the flap is vented if landing is the dominant maneuver and sealed if takeoff 
requirements govern the design. 

One advantage of the Kruger flap is that it can be applied to a wing with laminar 
flow technology. Since the Kruger flap is stowed along the lower surface of the wing, it 
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doesn’t disturb the flow over the critical upper surface of the wing in cruise flight (Fig. 
2.7). The presence of irregularities on the upper surface created by the trailing-edge of a 
retracted slat will trip a flow from laminar to turbulent and subsequently cause adverse 
effects on the performance. For this reason slats can not be used on laminar flow wings. 
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Figure 2.1 Lift Curves for Various High-Lift Airfoils 


Plain Flap Slotted Flap 



Figure 2.2 Trailing-Edge Flap Configurations 
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Variable Camber Kruger Linkage Fixed Hinge 

Figure 2.3 High-Lift Support Types 
(adapted from Ref. 7) 



Figure 2.4 Definition of Fowler Action 
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Figure 2.5a Usable Lift of Civil Transport Aircraft 
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Figure 2.5b Tail-Scrape Angle of a Typical Civil Transport Aircraft 
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Figure 2.6 Various Leading-Edge Devices 
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Figure 2.7 Surface Disturbances Created by Leading-Edge Devices In Cruising Flight 
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3.0 Performance Requirements 

The takeoff and landing performance of a subsonic civil transport aircraft is 
governed by the requirements listed in the Federal Airworthiness Regulations (FAR) Part 
25. These regulations specify minimum speeds, field lengths, and rates of climb that an 
airplane must maintain during a takeoff or landing maneuver. The high-lift system must 
ensure that the aircraft complies with these operating rules given basic fixed parameters 
such as wing loading and thrust loading. 

3.1 Takeoff 

Takeoff performance is characterized by the balanced field length and climb 
gradient of an aircraft. As illustrated in figure 3.1, the maneuver is divided into three 
segments: ground roll and rotation, lift-off and first segment climb, and second segment 
climb. There is also a final climb segment at 1 500 feet, but this manuever is performed 
with retracted high-lift devices. 

The ground run ends when the aircraft rotates at a specified speed and then lifts 
off after the speed reaches 1.1 (1.05 with one engine out) times the minimum unstick 
speed V MU , where V MU represents the minimum airspeed at which the airplane can safely 
lift off and continue the takeoff with its critical engine inoperative. This maneuver is 
generally a function of the maximum usable lift (not the absolute maximum lift), because 
an aircraft's rotation angle and, thus, angle of attack is restricted by the fuselage's tail- 
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scrape angle. The effects of this limitation must be taken into consideration when 
choosing a high-lift configuration, since the addition of another flap element may be 
required to increase the usable lift. 

The aircraft must now reach its takeoff climb speed V 2 before it arrives at the 
screen height of 35 feet. This speed must equal or exceed 1.1 times the minimum control 
speed and 1 .2 times the minimum dynamic stall speed. Generally the latter constraint 
corresponds to 1.13 of V S] , the stall speed in steady flight. This implies that the lift 
coefficient must be equal to or less than C Lmax j 1.13 2 . Thus, a high maximum lift 

coefficient is essential to obtain a low takeoff climb speed. 

The second segment climb begins when the landing gear is retracted. At this 
point the aircraft must maintain an airspeed greater them V 2 and a climb gradient greater 
than 2.4% for a twin-engine configuration, 2.7% for a tri-engine configuration, and 3.0% 
for a quad-engine configuration with one engine inoperative. This climb gradient is of 
great importance to a high-lift system, because it governs the necessary efficiency of the 
device. The climb angle is a function of the lift-to-drag ratio for an aircraft with a 
specified thrust loading: 


siny 


T_ J_ ' 
W~ L/D 


( 6 ) 


Therefore, an aircraft's climb performance can be improved with a high-lift system 
having a high lift-to-drag ratio. However, this presents a problem in the design of the 
overall system, for the first and second segments in climb require opposing capabilities. 

The high usable lift and maximum lift requirements of the first climb segment 
demand the use of a leading-edge device, more flap elements, and/or higher flap 
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deflection angles. This tends to degrade the climb performance needed during the second 
segment due to a decrease in the lift-to-drag ratio. Consequently, it is important to have 
the capability to accurately analyze the trade-offs between various high-lift configurations 
in order to optimize the aerodynamic performance of an aircraft so that it satisfies the 
FAR Part 25 requirements for both takeoff and landing. 

3.2 Landing 

The landing performance of a civil transport aircraft consists of an approach, 
transition from threshold to touchdown, and a braked ground run. The aircraft 
approaches at a glide slope angle of 3° and must maintain minimum approach and 
landing speeds as shown in figure 3.2. 

The FAR Part 25 requirements specify that the approach speed V A is at least 1 .3 
times the minimum dynamic speed V MS or 1.22 times the minimum stall speed V s in 
steady flight. Therefore, a low approach speed necessitates a high maximum lift 
coefficient; however, the pilot's ground visibility may limit the aircraft’s angle of attack 
during approach. This must be kept in mind during the design phase, because a leading- 
edge device may push the angle of attack at maximum lift, and consequently the angle of 
attack at 1 .22-V S[ , well beyond a safe pitch attitude for visibility. If this occurs, a more 
complex flap system may be required to shift the lift curve up in order to increase the 
usable lift coefficient. 

Climb gradient is also a factor during the landing maneuver if an aircraft must 
abort an approach. In general, a lower L/D is desired so that a higher thrust setting can be 
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used for improved handling and response. However, a climb gradient of 3.2% must be 
maintained with the flaps deployed, the gear down, and all engines operating if an aircraft 


has to abort its landing, and as a result a high lift-to-drag ratio is required. Once again 
trade-offs are necessary, and the high-lift system must be iterated upon to optimize the 
aircraft's performance. 
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4.0 Design Constraints 

The objective of a high-lift design is to produce the simplest configuration which 
meets all of the performance requirements. However, this is a difficult job because many 
limitations hinder the simplicity of the system making its optimization complex. 

The wing of a civil transport aircraft is designed to maximize M L/D for a given 
payload and mission. Subsequently, the wing area, sweep angle, aspect ratio, twist, and 
thickness are generally set in the early stages of a preliminary design. One goal of this 
research is to improve the interaction between the low speed and cruise aerodynamics 
during the early design phase so that high-lift performance may be considered when 
designing the wing. This may lead to improved low speed performance, yet the shape of 
the cruise wing will still impose many restrictions on the design of the high-lift system. 

The leading and trailing-edge devices may have constraints placed on their chord, 
span, or thickness ratio. Given these limitations the designer may have to increase the 
complexity of the system to comply with performance requirements. For example, if the 
sweep of the wing degrades the efficiency of a high-lift system to the point where a 
single-slotted flap is insufficient in producing the required lift, a second flap element may 
have to be added to increase the effective camber. 

The chord length of a slat or flap can be restricted by spar location or internal 
storage space within the wing. This can happen if a wing suffers from excessive twisting 
or bending and needs increased structural stiffness. Flap chord can also be limited if an 
aircraft requires a larger internal fuel volume. In either case, the front or rear spar 
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location may have to move leaving less room for a high-lift device. In addition, valuable 
space behind the rear spar may be restricted by landing gear, allowing little room for a 
retracted flap. 

The maximum thickness of a flap element is also a parameter which is often 
restricted for a designer, because most cruise airfoil shapes are very thin in the trailing- 
edge region. Although thickness effects of flaps are minimal (especially when employing 
a leading-edge device), excessively thin flaps can reduce C, ^ due to flow separation 
near the leading-edge of the flap. Slender flaps also have a relatively low structural 
stiffness. Consequently, their panel weight tends to be higher for a given stiffness, and 
they usually require more supports to maintain optimum gap settings. 

Many factors determine the span or continuity of a flap. For instance, inboard 
ailerons are sometimes necessary for high speed flight, or a thrust gate may be needed 
because of a short engine mount. This is a very important issue, because nothing is more 
detrimental to a flap's effectiveness than a discontinuity in its span. A trade-off study 
should be performed to analyze the effects of engine location on landing gear weight and 
flap efficiency, since the increase in landing gear weight may be far less than the increase 
in flap weight if an additional elerfient must be added to the trailing-edge to compensate 
for a flap cut-out. 

One last problem encountered in high-lift system design is that the shape of the 
leading-edge device is governed by the configuration type. The leading-edge shape of a 
slat is defined by the cruise wing. This leaves only the shape of the slat cove and the nose 
of the fixed leading-edge as a design variable. In contrast, a Kruger flap offers versatility 
in the contour of the leading-edge for optimum performance since it is stowed along the 
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lower surface of the wing when not in use. Because of this, some argue that a Kruger flap 
has the potential for a higher L/D and C L . 7 Nevertheless, a Kruger flap is a very 
complex device and is generally heavier than a slat. These are the kinds of decisions that 
are the target of this systematic evaluation algorithm. This research project is intended to 
develop tools that can be used to better understand the tradeoffs quantitatively so that the 
optimal configuration may be selected. 
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5.0 Computational Modeling 


Today it is virtually impossible to conduct an extensive experimental program to 
generate a high-lift database as was done in the United Kingdom during the 1970's and 
80's. 3 * 5 * * However, advances in computational fluid dynamics have made it possible to 
accurately, yet inexpensively, predict the flow field around multi-element airfoils. 
Therefore CFD was used as the primary tool in this research project along with some 
experimental data. The following three viscous-flow solvers for multi-element airfoils 
were used to generate an aerodynamic database: 

1 . INS2D 8 - 9 This code solves the incompressible Navier-Stokes equations on 

structured overset meshes. The code was developed by S. Rogers at the 
NASA Ames Research Center and is widely used in industry for high- 
lift research. 10 * * 

2. NSU2D 1112 This code solves the compressible Navier-Stokes equations on 

unstructured meshes. The code was developed by D. Mavriplis at 
ICASE and is also widely used in industry. 13 

3. MSES 14 - 15 This code solves the Euler equations and integral boundary-layer 

equations simultaneously. It was developed by M. Drela of MIT and at 

the high-lift workshop CFD challenge it was shown to be the most 

accurate viscous/inviscid interaction method. 


* 

These three codes are fairly well validated and have been shown to predict 

changes in the forces and moments due to small differences in a flap configuration. 16 The 

aerodynamic data obtained with the codes was augmented with experimental data 

whenever available. The advantage of this hybrid (i.e., computational as well as 

experimental) approach is that innovative concepts can be analyzed at flight Reynolds 
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numbers and included in the database much more quickly and inexpensively using CFD 
than using wind-tunnel experimentation. 

Other instances required the use of experimental data because of current CFD 
limitations. For example, present computational analysis of Kruger flaps is not possible, 
because the separated flow region behind the flap cannot be modeled accurately. As a 
result, wind-tunnel data was used to generate the Kruger flap portion of the database. 

Special note should also be made that while CFD is able to accurately compute 
the lift of relatively complex high-lift airfoils, the prediction of drag for these airfoils still 
remains a challenge. Some engineers argue that the errors involved in predicting drag 
arise from the inadequacies of the turbulence models used by Navier-Stokes codes. 
However, this may be only partially true as explained by Vinh et al . 17 They show that 
even the prediction of drag for multi-element airfoils with attached flow can be 
inaccurate. They suggest that these errors arise from the method of integration used to 
calculate the drag force and demonstrate that wake integration techniques provide 
improved results over surface integration. 
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6.0 Aerodynamic Database 

Computational Fluid Dynamics was used to construct a database consisting of the 

NLR-7301 two-element and the Douglas LB-546 three-element airfoils. However, the 

new design methodology is easily extensible, and new systems can be added to the 

database without difficulty. The two-element and three-element airfoils shown in figure 

6. 1 were considered because they currently offer high aerodynamic performance. They 

have also been tested extensively in wind-tunnels and are therefore well documented. 

INS2D was used principally to obtain CFD results of the NLR-7301 two-element 

airfoil. A sample grid of the airfoil used for the calculations is illustrated in figure 6.2. 

Figure 6.3 shows excellent agreement between INS2D’s results and the experimental 

data 18 for a pressure distribution at an angle of attack of 6° and Reynolds number of 

2.5 lxl 0 6 . The lift curve for the two-element airfoil can be seen in figure 6.4. The data 

compares well at low to moderate angles of attack; however, a slight deviation exists at 

higher angles and at maximum lift. The computational results exhibit an extended stall 

that continues a few degrees past the wind-tunnel data. This apparent overprediction of 

* 

C. has also been encountered by Lin and Dominik. 19 INS2D seems to consistently 

‘max 

overestimate the stall angle resulting in a higher value for maximum lift. 

Figure 6.5 contains the computed and experimental drag polar of the NLR-7301 
airfoil, and again the results show excellent agreement between INS2D and the wind- 
tunnel data. Also shown in the figure is a drag polar which was fit to the computational 
data using the RAF parabolic polar estimation method developed by Lean and Fiddes. 20 
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The RAF method was used to construct a drag polar for each high-lift configuration in 
order to develop the drag prediction method outlined in section 7.1 . 

INS2D was also used to obtain data of the Douglas three-element airfoil, and an 
example of a grid used in the computations is illustrated in figure 6.6. Two predicted 
pressure distributions of the airfoil at an angle of attack of 8.1° and Reynolds number of 
9xl0 6 are shown in figure 6.7. The difference in the two cases reflects the modeling of 
the farfield conditions, since grids of the airfoil with and without the wind-tunnel wall 
were generated to compare with the experimental data. As reported by Kusunose and 
Cao 21 , including a wind-tunnel wall improved the agreement with the experiment. This is 
especially evident over the slat, where the case with a wind-tunnel wall did much better in 
predicting the pressure. 

Figures 6.8 and 6.9 illustrate the lift curve and drag polar respectively of the 
Douglas airfoil. The computational results in these plots represent only the airfoil 
without wind-tunnel walls, so the results are slightly skewed from the experimental data. 
Regardless of this discrepancy, the data was used in the construction of the database for 
this configuration since it more accurately models an airfoil in steady flight. 
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Figure 6.1 Two and Three-Element Airfoil Geometries 



Figure 6.2 Sample CHIMERA Grid of the NLR-7301 Two-Element Airfoil 




2 1 

Figure 6.3 NLR-7301 Pressure Distribution (a=6°, Re=2.51xl0 6 ) 



a (deg) 


Figure 6.4 NLR-7301 Lift Curve (Re=2.51xl0 6 ) 
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7.0 Aerodynamic Module 

The newly developed high-lift module predicts the aerodynamic performance for 
a given high-lift system to be used within a takeoff and landing optimization routine such 
as ACSYNT. It calculates lift, drag, pitching moment, and maximum lift for a given 
high-lift system and flight conditions. 

The high-lift module first calculates two-dimensional aerodynamic characteristics 
from equations based on CFD and wind-tunnel data. These aerodynamic characteristics 
consist of sectional lift curve, profile drag polar, and pitching moment which are 
calculated from basic airfoil geometry. A modified lifting-line method based on 
Weissinger's theory 22 is then used to determine the total aerodynamic coefficients of the 
wing. Thus total lift, drag, and pitching moment can be computed faster than with a 
panel method while incorporating (based on strip theory calculations) the viscous effects 
of slotted high-lift devices. 

7.1 Two-Dimensional Aerodynamics 

■% 

Theoretical and empirical techniques were applied to develop equations based on 
the high-lift configuration type, c, t/c, c LE/XE , 8 LE/rE , a, a 0 , Q a > s ” anc * s * ^ rom two " 
dimensional database in order to calculate a n , C, , C d , and C m . 

u a profile o 
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7.1.1 Lift 

A flap is used to produce added lift by increasing the effective camber of an 
airfoil. Thus, the total lift of an airfoil with extended high-lift devices is represented by 
the lift of the clean airfoil plus the incremental lift created by the leading and trailing- 


edge flaps: 


C '= C ' c | OT 7 + AC 'TE +AC 'LE 


(7) 


The lift coefficient for the clean airfoil at a given angle of attack is provided by 
ACSYNT. The increment in lift coefficient is a function of the flap configuration, 
geometry, and deflection angle. The additional lift produced by a single-slotted trailing- 
edge flap is shown below as defined in DATCOM: 23 


AC 1t =C, <x s 5 te — 

*TE “clean 6 111 c 


( 8 ) 


The lift-effectiveness, a 5 , is determined for a specified flap configuration type with the 
general equations shown below: 


<x s = — 


-TE 


+ X 


a S theory ^t/c^visc + Q 


TE 


a 


^ theory 


n 




C TE f +s j n -l I^IE. 


X = [l + a, • tan -1 (a 2 -5 TE x i )] fcf 


f Cf =l-b, 


"TE 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 
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f,/ c — 1 + c t • t/c (13) 

fvisc=d, (14) 

where a,, a 2 , b,, c„ d,, and x, are empirical coefficients and were determined to be as 
follows for both the NLR and Douglas flaps: 

a, = 0.2, a 2 = -5.2, b, = 0.55, c, = 1.04, d, = 0.785, x, = 5 
The expression for the theoretical lift-effectiveness, (* 5 Iheoiy , is derived for a bent 
flat plate in an inviscid flow and is a good approximation for small flap deflections. 24 It is 
a function of the flap chord, c re , only, but an empirical airfoil thickness factor, f l/c , and 
viscous factor, f vjsc , have been added to account for Reynolds number by scaling a 5( 
with CFD or experimental data. The separation factor, x, is used to reduce the lift- 
effectiveness, a 5 , at higher flap deflection angles. The empirical factor is roughly unity at 
small flap deflections but decreases at larger deflections to account for strong viscous 
effects and flow separation over the flap. This factor is dependent on the configuration 
type and the Reynolds number of the CFD or experimental data. Examples of the 
separation factor and lift-effectiveness curves for the NLR-7301 airfoil are illustrated in 
figures 7.1 and 7.2 respectively. 

It should be noted here that the lift increments of the NLR flap and Douglas flap 
calculated with this method are both higher than the lift increment of a double-slotted flap 
calculated via DATCOM. 23 This illustrates the importance of the present methodology, 
since it is obvious that the performance of current high-lift systems can no longer be 
predicted with methods that are based on outdated technology and low Reynolds number 


data. 
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The primary task of a slat or Kruger flap is to extend the lift curve by increasing 
the stall angle and as a result increase a wing's maximum lift capability. However, a 
leading-edge device will slightly alter the lift of a wing. This lift increment is often 
negligible and may be disregarded in the preliminary design phase, but if the effect is 
desired, linear theory may be used for a good first approximation. The additional lift 
produced by a leading-edge device is calculated from: 


c' 

AC| LE _ ^'adean ^theory ^ LE c 


(15) 


where Roshko 25 shows that the theoretical lift-effectiveness of a leading-edge device is: 



The extended chord length, c' , is calculated from geometric properties of the 
leading and trailing-edge devices. This increase in chord is produced by so-called Fowler 
action and is a very important characteristic of a high-lift device. It can dramatically 
improve the lift performance of an aircraft without producing excess drag. Examples of 
Fowler action are shown below for an airfoil with a leading-edge device and a single or 


double slotted flap. 

* 

Single slotted flap: 

c' =c + f-s t +d-c LE (17) 

The Fowler action functions, f, of the flap and, d, of the slat can be specified by the user, 
or the default functions may be used. The flap function is based on the assumption that 
full Fowler motion is reached at a deflection angle of 45° and is specified as: 
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f = 


fsin(28 XE ) for 5 TE <45° 


1 


for 5 te > 45° 


(18) 


The leading-edge extension function assumes an increase in Fowler action of 70% of the 
leading-edge chord for a takeoff position of § LE = 20° and 90% for a landing position of 


5 le =30° 


d = sin(25 LE 09 ) (19) 

Generally a sealed leading-edge device is employed during takeoff and a vented device is 
used when landing. 

Double slotted flap: 

c' = C + f[ -s, + f 2 ■ s 2 +d -c LE (20) 

The shroud length of the fore flap element is labeled s 2 , and the functions f, and f 2 denote 
the Fowler motion of the vane and the aft flap respectively. These Fowler motion 
functions shown below specify full fore flap deployment at 45°, beginning deployment of 
aft flap at 20°, and full aft flap deployment at 50°: 


H 

sin(25 TE ) 

1 

for 

for 

5 te 

8 te 

<45° 

>45° 

(21) 


0 


for 

5 te < 20° 


f 2 =- 

sin[3(5 XE 

-20)] 

for 

20° < 8 xe <50° 

(22) 


1 


for 

8 te > 50° 



The lift curve of the high-lift airfoil can now be determined using the above 
information. The local lift curve slope for the airfoil is given by: 
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C, = X C, 

‘a K ‘a 


clean c 


( 23 ) 


Note that the separation factor % appears once again; this time it accounts for a deviation 
in the lift curve slope from 2n rad' 1 . Finally, the angle of attack at zero lift for a flapped 
airfoil can be calculated from: 

C, 

(24) 


a » = -c 


where the local lift coefficient at zero angle of attack is determined as shown below: 


C, = C, -+ac 1tf + ac Iif 

■o °clean c TE LE 


(25) 


7.1.2 Drag 


A database of minimum drag coefficient, Ch , lift coefficient at minimum drag, 
C, , and profile drag constant, k_, was constructed using the RAF parabolic polar 

min r v 

estimation method. Once this database was established, equations for the increments in 
C H , C, , and were determined by fitting curves to the data. The following equations 
represent the increments associated with a single-slotted trailing-edge flap: 


AC d . =a,5 TE !< 


C TE 1 2 c 

V 0.3 


(26) 


AC lmin TF “( b l 5 TE ~ b 2 5 TE X3 
TE 


C TE 1 4 c 

0.3 


(27) 


A ^ p TE _ ( C l b TE +C 2 b TE 2 )l 


Ct E l ^ c 
03 


(28) 
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where the empirical constants and exponents given below are for both the NLR and 
Douglas flaps: 

a,=0.038, b,=2.3, b 2 =2, c,=0.00012, c 2 =0.0097. 


x,=1.74, x 2 =1.4, x 3 =4, x 4 =-0.17, x 5 =2.55 
The incremental drag coefficient for a leading-edge flap is extremely dependent 
on the configuration type, deflection angle, and gap setting. Therefore, no general 
empirical method for calculating drag increment has been developed. Instead CFD data 
is used explicitly to give drag increment based on the configuration tested. The Douglas 
slat is used here to give drag polar information for both the takeoff and landing settings: 
Takeoff (5 le =20°): 




= 0.0013— 


‘LE 


c' 

AC, . = 0.46— 

*min LE C 


(29) 

(30) 


Ak D = 0.00772 -7 

f LE C 

Landing (8 LE =30°): 


(31) 


AC d . =0.0074 — 

dmm LE C 


c' 

AC, . =0.76 — 

min LE C 


Ak D =0.00381- 

P LE C f 


(32) 

(33) 

(34) 


The drag polar components of the high-lift airfoil are calculated from the 
contributions of both the leading and trailing-edge devices and the clean airfoil: 
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Ch . = c. 


'd min d min c|ean c 


+ AC d . + AC d . 

a min^ amm LE 


Ci . = C | . + AC | . + AC | 

mln mm clean C mm TE 


k p = k W + AI %. +Ak ^ < 37 

These coefficients are then used to construct the following airfoil drag polar which is 
used to determine profile drag at each section: 

C d ,, =C d . +k (C, -C, . ) 2 (38 

^profile “mm P v 1 •min 7 v 


7.1.3 Pitching Moment 


The pitching moment about the aerodynamic center (or at zero lift) of the high-lift 
airfoil is calculated from the pitching moment of the clean airfoil and the contributions 
from the leading and trailing-edge devices: 


C = C 


m ° m °clean V C 


I +AC m +AC m 

C / °TE °LE 


The increment in the pitching moment coefficient at the aerodynamic center due 
to the trailing-edge device is: 


1 X cp c' 

AC m =AC, — 7 — 

°te ‘te14 c' c 


The position of the center of pressure, — — , for an airfoil with a single-slotted flap is 

c 


determined empirically as: 
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— ~ = 0.5 — a,5 TE x i (41) 

0 

The empirical coefficients a, and x, were determined to be as follows for both the NLR 
and Douglas flaps: 

a, =0.134, x, = 0.46 


The pitching moment increment at the aerodynamic center due to the leading-edge 
device is: 


AC mo “Cm. ^ LE I 

°LE 8 LE V C 


»\2 


(42) 


where the pitching moment effectiveness about the aerodynamic center is found from thin 
airfoil theory: 26 


C m =-^ sin 0 LE(l-COSe LE ) 


(43) 


and 0 LE is specified as: 


n -1 , ~ v LE 

0 LE = cos 1 - 2 


(44) 


A re f . 

The sectional pitching moment about an arbitrary reference point — — is 
measured from the leading-edge of the clean airfoil is determined as: 

C m Xref =C m 0 +(C l cosa + C d sina)(^-- 


A 


ac 

C J 


( 45 ) 


where C d includes both the profile drag and induced drag at the section. 
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The aerodynamic center of the high-lift airfoil is only a function of Fowler action, 
and its position relative to the leading-edge of the clean airfoil is found by determining 

the individual shifts in — — due to the leading and trailing-edge devices: 

c 


x ac 

c 



(46) 


where 
edge device. 



is the extended chord resulting from only the Fowler action of the leading- 


7.2 Three-Dimensional Aerodynamics 

A modified lifting-line method based on Weissinger's theory is applied to 
determine the spanwise loading across the wing. From this procedure a load distribution 
can be calculated for a wing of arbitrary planform from sectional lift curve slopes and 
angles of zero lift. Thus, the total aerodynamic coefficients of the wing can be calculated 
from a combination of two- and three-dimensional methods. This results in a prediction 
method which is relatively fast and reflects the viscous nature of slotted flaps. 

The Weissinger method accounts for partial span flaps/slats, aspect ratio, and 
Mach number in order to calculate lift, downwash angle, induced drag, and pitching 
moment. The geometry of the wing may be asymmetric about the centerline and can 
include sweep, taper, and twist. The method models the wing as a plate of zero thickness, 
but the planform and twist remain identical to the actual wing. Thus, partial span flaps 
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and Fowler action may be modeled as discontinuities in twist and an increase in chord 
length respectively. 

The chordwise load distribution at each span station is concentrated into a lifting 
line located at the wing's quarter-chord line. The theory requires that the quarter chord 
line is straight, yet a discontinuity is allowed at the plane of symmetry so that a swept 
wing may be considered. 

The method assumes that the lifting line and its trailing vortex sheet are 
continuous. However, discrete values for the circulation strength are determined only at 
the span stations along the quarter-chord line corresponding with control points. The 
number of control points is specified, and each point is placed at the three-quarter-chord 
line. Mathematically this implies that the lift curve slope is 27i rad' 1 , so a correction 
method must be incorporated for a lift curve slope that varies from the theoretical value. 
For further details please refer to the outline of this correction method in Appendix A. 

The general input parameters required by the method are a, AR, r| LE/rE , A, X, M, 
and S, and the following parameters must be specified at each control point: a 0 , C, , c, s vv . 
The chord is defined to be parallel to the freestream direction in the Weissinger method 
as shown in figure 7.3, so caution'should be taken to ensure that all aerodynamic data 
represent streamwise values. Simple sweep theory dictates that two-dimensional data 
should be modified from normal (n) to streamwise values as follows: 


z/c = (z/c) n -cos A 

(47) 

M = M n /cosA 

(48) 
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Re = Re n /cos 2 A 

(49) 

C, =(c, ) ’Cos A 

(50) 

'a V 'a/ n 

c m =(c m )„ COSA 

(51) 


Also, recall that a„ and C, were determined previously from sectional CFD and 
experimental data and are based on the appropriate Reynolds number. 

7.2.1 Lift 


The Weissinger method determines the local lift coefficient, C„ at each control 
point along the wing. From this load distribution the total lift coefficient of the wing, C L , 
is determined. The method also returns values for the total lift coefficient of a wing at 
zero angle of attack, C Lq , the lift curve slope of wing, C L(j , and the downwash angle 
distribution, s. 

Two test cases were used to validate the results of the Weissinger method. Figure 
7.4 shows one of these test cases. It is a part-span-flap wing/body configuration with 
variable twist, and figures 7.5 and' 7.6 compare the experimental data 27 of this model with 
predicted load distributions. The other test case is a plain swept wing and is shown in 
figure 7.7. Similarly, figures 7.8 and 7.9 compare the experimental data 28 of this 
configuration with calculated lift distributions. Both validations compare relatively well 
with the experimental data; however, there is a slight disagreement at the root of the 
wing. This discrepancy exists because a body was not modeled computationally in either 
test case, while both of the experimental models were wing/body configurations. 
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7.2.2 Drag 

Total drag is calculated from both profile drag and induced drag contributions. 
The Weissinger method determines only the induced drag coefficient, C D ,, so further 
calculations are required to determine the total drag coefficient of the wing. 

Two-dimensional drag polar data is used to calculate the total profile drag by 
integrating local profile drag across the wing: 


C v» = C an,m + M C ' atSta< ' 0n 1 

(52) 

1 y body j b/2 

C Dprofile = g J C dCdy+ g 

^ y body 

(53) 


where y body represents the y-coordinate of the exposed wing root. For a symmetric wing, 
this reduces to: 

2 b/2 

^Dprofjie = g J[Cd ml „+kp(C|-C lmlll ) 2 ]aJy (54) 

y body 

The total drag coefficient can now be found by summing the profile drag and induced 
drag coefficients: 

profile 


D induced 


( 55 ) 
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7.2.3 Pitching Moment 

The total pitching moment about a reference line which is normal to the 
freestream velocity is calculated as: 

b/2 

l c ”«, t , cdy 

’ body 

where the local pitching moment is calculated as shown earlier: 

C mxref =C mo +(C,cosa + C d sina)[^--^] 

Here, the position — — corresponds to a point on the reference line as shown in 
c 

7.10. 

7.3 Maximum Lift Prediction 

Maximum lift and stall angle are also calculated using both two- and three- 

dimensional methods. The maximum local lift coefficient at each spanwise station along 

* 

the wing is calculated from empirical equations generated from the two-dimensional 
database. Weissinger's theory is then used to determine the onset of stall by finding the 
critical section. This method is used to predict maximum lift in the preliminary design 
process, but other more sophisticated methods should be used in subsequent design 
phases. 


(56) 


(57) 


figure 



One suggested method that can be implemented into the present methodology is 
the pressure difference rule developed by Valerezo and Chin. 29 This method predicts the 
maximum lift of an airfoil based on the difference in the peak pressure coefficient and the 
pressure coefficient at the trailing-edge of an airfoil. The pressure difference rule 
compares chordwise pressure distributions scaled by their corresponding peak pressure 
coefficients with a prescribed pressure difference at maximum lift to determine the onset 
of wing stall for a specified Mach and Reynolds number. This correlation has shown 
remarkable results and has demonstrated that the method works well for both clean and 
multi-element airfoils alike, assuming the leading element stalls first. 

7.3.1 Two-Dimensional 


Empirical equations based on CFD and wind-tunnel data are used to calculate 
maximum lift at each specified wing station. The general input parameters used to 
calculate maximum lift are: high-lift configuration type, c, t/c, c LE/TE , s,, s 2 , 5 LE/TE , r| LE/TE , 
A, M, and Re. 

The maximum local lift cofefficient is calculated from the maximum lift of the 
clean airfoil and the maximum lift contributions from the leading and trailing-edge 
devices: 


C, = C, 


c' 


+ AC i + AC i 


max 'max c | ean c 'max^E 'maxjR 


( 58 ) 


- clean maximum lift coefficient 


clean 
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AC, 


AC, 


X LE 


TE 


- lift coefficient increment due to the leading-edge device 

- lift coefficient increment due to the trailing-edge device 


The increment in maximum lift from a leading-edge device results from a shift in the stall 
angle of the airfoil: 


AC, = C, Aa max 

‘maxLE ‘a niaA 


= 3 1 + a 2 


'LE 


+ a- 


'LE 


+ 


"LE 


(59) 

(60) 


The shift in stall angle is only a function of configuration type and leading-edge chord 
length since the maximum lift performance is fairly insensitive to the angle of the 
leading-edge device if the gap and overhang are properly optimized. 5 30 As mentioned 
earlier, this shift in angle is constant regardless if the trailing-edge flaps are stowed or 
deployed. Figure 7.1 1 illustrates Aa max , and table 7.1 gives values for the empirical 
constants a„ a 2 , a 3 , and a 4 . 

The increase in maximum lift created by trailing-edge flaps is determined using a 


type of maximum lift-effectiveness parameter, C, : 

^inax 


AC 


— Ci r 5 tp 

c 


m ax T c lx U TE 

max TE °max 


C, 8 max "( a > +a 2 5 TE X ‘) f cf 


, C TF ' X2 

f = — 

Cf ' 0.3. 


(61) 

(62) 

(63) 


where the empirical constants used to calculate the increment in maximum lift for single- 


slotted flaps are: 
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a,=l 53.4, a 2 =-l 51.8, x,=0.018, x 2 =0.16 
7.3.2 Three-Dimensional 

Weissinger's theory is applied to determine the lift distribution across the wing, 
and the maximum lift is predicted using a critical section approach. It is assumed that the 
onset of stall occurs when one wing section reaches its maximum lift as illustrated in 
figure 7. 1 2, yet the wing continues to lift past this critical angle of attack as reported by 
Murillo and McMasters. 11 The maximum lift of the wing is evaluated as: 

C Lm „ =C L (a crit ) + f(Re) (64) 

where C L (a cnt ) corresponds to the lift coefficient at the critical angle of attack, and the 
function f(Re) is Reynolds number dependent and is assumed to be 10% of C L (a crit ) at the 
present. 

The polynomial method for optimization 32 is used to find the critical angle of 
attack since sectional lift curves may be non-linear. The difference between the local lift 
coefficient and the local maximum lift coefficient is found at each spanwise station from 

the provided load distribution. The minimum difference, (C. - C, , is determined, 

V max ) m j n 

and the objective function is defined as: 

The polynomial optimization method is then used to minimize this function in order to 
determine the critical angle of attack. 




50 



Sealed Slat 

Vented Slat 

Sealed Kruger 

Vented Kruger 

a, 

0.0068 

0.017 

0.0068 

0.017 

a 2 

0.81 

1.28 

0.81 

1.28 

a 3 

-1.88 

-1.35 

- 0.37 

0.32 

a 4 

-0.38 

-5.5 

-0.357 

-5.3 


Table 7.2 Variables Used to Calculate the Maximum Lift-Effectiveness 

of Leading-Edge Devices 
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Figure 7.1 NLR-7301 Separation Factor 
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Figure 7.2 NLR-7301 Lift-Effectiveness 
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2-D 

Reynolds Averaged N-S Method 


3-D 

Modified Lifting-Line Method 


<k 



Section A-A 
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▼ 


i A 


Figure 7.3 Overview of Aerodynamic Analysis 



53 


Wind-tunnel measurements using a variable twist wing (VTW) model in the NASA 
Langley 14- by 22-Foot Tunnel at a Reynolds number of lxlO 6 and lift coefficient of 0.6. 



* Cylindrical body length 

Figure 7.4 Part-Span-Flap Deflection Test Case (from Ref. 27) 



y/(b/2) 


Figure 7.5 Part-Span-Flap Spanwise Load Distribution 
VTW7S 1 Configuration, a= 1 1 .4° 
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Figure 7.10 Procedure for Calculating Pitching Moment 
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Figure 7.11 Increment In Stall Angle Due to Various L. E. Devices 
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Spanwise Station, r| 

Figure 7.12 Maximum Lift Estimation Using Critical Section Approach 




58 


8.0 Weight Module 

The weight of the high-lift system is an important factor in the preliminary design 
phase of an aircraft. The system weight is very much governed by the aerodynamic 
loads, the structural materials, and the structural stiffness requirements. Other important 
factors that may affect the overall weight of the high-lift system are fail-safe design 
requirements, system complexity, and system reliability. However, it is important to note 
that information on the weight and the aerodynamic performance of a high-lift device is 
not sufficient to determine the optimal system. The optimal system depends on the 
effects of the high-lift system on the weight and cost of the complete aircraft. For 
instance, a variable camber Kruger flap may be heavier than a slat; yet, as a result of a 
slightly higher maximum lift, it may produce a substantial increase in maximum payload 
for a given aircraft geometry and mission profile. 

This estimation method was developed to offer a parametric weight sensitivity 
analysis for weight and cost trade-off studies. It uses empirical equations based on 

historic data to calculate incremental weights due to leading and trailing-edge devices. 

% 

The method predicts leading and trailing-edge surface, support, support fairing, actuation, 
and fixed structure weights for a given high-lift system. As a result, changes in the 
technology of materials, geometry, mechanisms, or actuation can be reflected through the 
scaling of the various weight components. 

The equations are based solely on stowed flap area, so they require minimal 
geometric data to evaluate system weight in the preliminary design phase. This 
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represents a reduction in complexity from existing weight estimation methods such as 
those by Anderson et al. B and Oman. 34 These methods, which are based on numerous 
high-lift parameters, are very difficult to correlate with the limited data that is available, 
and they often result in trends that are not realistic. Thus, a simpler approach was 
adopted to construct a fast, robust, and fairly accurate weight prediction method for 
multi-element wings. 

8.1 Trailing-Edge Flap Weight 

Flap panel, support structure, support fairing, and actuation weight correlations 
are shown below. These weights are based on stowed flap area for a given flap 
configuration and support type. The specific weights a„ a 2 , a 3 , and a 4 are given in tables 
8. 1 - 8.3. Trailing-edge weight correlations are shown for several aircraft in figure 8. 1 . 

Surface Structure: 

^panel“ a r^TE (^6) 


Support Structure: 


^support a 2 ’ ffow 


5 TE 


f 


fow 


0.47 + 0.53 


^ fow 

0.5 


(67) 


for linkage and external hinge supports 


1 


for track supports 


( 68 ) 
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Rro» = m 

b TE 

where f fow is a function of Fowler motion and is used to scale the weights of linkage and 
externally hinged supports and fairings. The term is the ratio of Fowler area to the 
area of the trailing-edge flap. 

Support Fairing: 

Wfairing = a 3 ' ffow ' ^TE (70) 

Actuation: 

W acUl ation = a 4 ‘ ^TE (71) 

Total Trailing-Edge Flap Weight: 

Wjg — w panc , + W SU pp 0rt + Wf airin g + W actuation (72) 

The fixed trailing-edge weight of the wing is also determined in order to calculate 
the total wing weight. This weight is estimated using the equation developed by 
Anderson et al.: 35 

W flxedTE = 1.34(s fixedTE ) > 2.0-S flX£dTE (73) 

where S flxedTE is the area of the fixed trailing-edge and is calculated as: 

^ fixed T E = S TE gross “ ( S TE + S ail + 05 ' S sp) ( 74 ) 
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The gross trailing-edge area of the wing is represented by , and the aileron and 
spoiler areas are denoted as S ai , and S sp , respectively. 

8.2 Leading-Edge Flap Weight 

As with the trailing-edge flaps, the leading-edge device component weights are 
based on nested flap area as shown below. The specific weights a 5 , a 6 , and a 7 are given in 
table 8.4, and correlations for leading-edge weights of several airplanes are plotted in 
figure 8.2. 

Moving Panels: 

W m ™„ s =a 5 S LE (75) 

Actuation: 

^actuation = a 6 '^LE (76) 

Fixed Leading-Edge: 

Wfixed = a 7 'S L e (77) 

Total Leading-Edge Flap Weight: 

W LE = w„ 0ting + W, ao>lion + W 0>t( , (78) 
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The above methods for predicting leading-edge device, trailing-edge device, and 
fixed trailing-edge weight along with the methods for estimating aileron and spoiler 
weight (shown in Appendix B) will be used in conjunction with the wing box weight 
estimation method developed by Ardema et al. 36 to determine the total weight of the wing 
within ACSYNT. 

This prediction method is intended to be used only in the preliminary design 
phase of an aircraft. However, the methodology lends itself to an improved weight 
estimation method for subsequent design phases, since experimental data shows 37 that 
flap loading is independent of angle of attack and is only a function of flap deflection. 

The proposed idea involves changing the routine from a 0 ,h order to a 1 51 order predictive 
method so that CFD results can be used in conjunction with slender beam analysis to 
determine the weight of a configuration. Thus, the chordwise load distribution supplied 
by CFD and the spanwise load distribution provided by the Weissinger method can be 
used to determine the weight of a high-lift system for a given flap placard speed 
(maximum speed allowed with flaps deployed) with a specified support type and number 
of supports. 

Once this method is developed, the loading information provided by CFD can also 
be used as an input in the calculation of wing box weight. This will provide valuable 
information for calculating the torsional loads produced by the flaps. Consequently, the 
structural analysis of the wing box will include information regarding the flap 


configuration and Fowler action. 
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Specific 

Weight 

Single-Slotted 

Fixed 

Vane/Main 

Articulating 

Vane/Main 

Main/Aft 

Triple- Slotted 

Flap Panel (a,) 

2.7 

3 

3.5 

4.8 

5.5 

Support (a 2 ) 

3.0 

3.2 

3.8 

4.7 

5.6 

Fairing (a 3 ) 

1.0 

1.0 

1.15 

1.3 

1.4 

Actuation (a 4 ) 

2.2 

2.2 

2.3 

2.4 

2.5 

Total 

8.9 

9.4 

10.75 

13.2 

15 


Table 8.1 Specific Weights of Trailing-Edge Flaps with Hooked Track Supports 

(lb / ft 2 of stowed flap area) 


Specific 

Weight 

Single-Slotted 

Fixed 

Vane/Main 

Articulating 

Vane/Main 

Main/Aft 

Triple-Slotted 

Flap Panel (a,) 

2.7 

3.0 

3.5 

4.8 

5.5 

Support (a 2 ) 

1.5 

1.6 

1.9 

2.4 

2.8 

Fairing (a 3 ) 

0.11 

0.11 

0.13 

0.14 

0.16 

Actuation (a 4 ) 

2.0 

2.0 

2.1 

2.2 

2.3 

Total 

6.31 

6.71 

7.63 

9.54 

10.76 


Table 8.2 Specific Weights of Trailing-Edge Flaps with Link/Track Supports 

(lb / ft 2 of stowed flap area) 


Specific 

Weight 

Single-Slotted 

Fixed 

Vane/Main 

Articulating 

Vane/Main 

Main/Aft 

Flap Panel (a,) 

2.7 

3.0 

3.5 

4.8 

Support (a 2 ) 

1.1 

1.2 

1.8 

2.2 

Fairing (a,) 

0.28 

0.28 

0.29 

0.29 

Actuation (a 4 ) 

0.9 

1.0 

1.3 

1.5 

Total 

4.98 

5.48 

6.89 

8.79 


Table 8.3 Specific Weights of Trailing-Edge Flaps with External Hinge Supports 

(lb / ft 2 of stowed flap area) 
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Flap Type 

Fixed Camber 
Kruger 

Variable Camber 
Kruger 

Slat without 
Slave Tracks 

Slat with 
Slave Tracks 

Movable L. E. (a 5 ) 

4 

5.1 

7.7 

8.5 

Actuation (a 6 ) 

2.8 

2.8 

2.1 

2.1 

Fixed L. E. (a 7 ) 

6 

8.5 

4.5 

4.7 

Total 

12.8 

16.4 

14.3 

15.3 


Table 8.4 Specific Weights of Leading-Edge Devices 
(lb / ft 2 of stowed flap area) 
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Figure 8.1 Trailing-Edge Weight Correlations for Various Aircraft 
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Figure 8.2 Leading-Edge Weight Correlations for Various Aircraft 
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9.0 Cost Module 

The impact that a high-lift system has on the total life cycle cost of a transport 
aircraft is a complex function of the costs to develop, build, purchase, operate, maintain 
and dispose of the airplane. As described by Gelhausen 38 , the total cost of the system is 
best captured in the direct operating cost (DOC) where the costs to own and operate the 
aircraft are combined. The DOC includes fuel, airframe, crew, and maintenance expenses 
and is the most important parameter used when comparing civil transport aircraft with 
equivalent missions. A full mission and market analysis is required to show the effects of 
improved high-lift system performance. This includes improved fuel economy, the 
increase in range and payload, and the increase in units sold. However, the cost to 
manufacture and maintain the system must be analyzed in order to understand the 
complete picture. For instance, a system that is more complex but generates a higher 
maximum lift may result in improved cruise efficiency and reduced fuel expenditures, but 
at the same time may cause increased airframe and maintenance cost. It is important to 
not only provide the aerodynamic and weight information, but also the cost information 
to the synthesis process. This methodology incorporates system cost so that careful 
tradeoff studies can be conducted to determine the optimal high-lift system for a new 
aircraft. 

The system cost is governed by the weight, W, and the part count, PC of the high- 
lift device : 39 


Cost = a r W-PC x ' 


( 79 ) 
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[ 1.888 1 trailing - edge device 

a ' [1.7339 leading - edge device 

x,=0.7 

Here the constant, a, accounts for hourly labor costs and type of material. Historically, 
weight has been the most important parameter, but recently more information regarding 
part count has also become available. For instance, Boeing simplified the flap system on 
the new 737 from a triple-slotted to a double-slotted design. The new flap system has 
37% fewer parts and 33% fewer bearings. 40 This reduction in part affects cost (note, 
however, that the reduction in cost is reported to be significantly less than 30% indicating 
a value for the part count exponent, x„ significantly less than unity) and, consequently, 
the newly developed cost module incorporates this factor as well as system weight. 

The total part count of a system is found by summing up the various components 
of the device. The part count of a trailing-edge flap consists of the flap panel, support, 
support fairing, and actuation: 

PC TE = PC panel + PC SU pport + Pairing + PC actuation (80) 

The part count of a leading-edge device is the sum of the movable leading-edge, 
actuation, and fixed leading-edge:* 

PC LE = PC mova bie + PC actU ation PC fixed 

Each category is a function of the flap configuration and support type only. The total part 
count is used as a measure of overall system “complexity”, while the size of the device 
and complexity associated with Fowler action are reflected in the weight of the flap. 
Values for part count components of trailing-edge devices are given in tables 9.1 - 9.3, 
and component values for leading-edge devices are given in table 9.4. 
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Part Count 
Component 

Single-Slotted 

Fixed 

Vane/Main 

Articulating 

Vane/Main 

Main/Aft 

Triple-Slotted 

Flap Panel 

600 

750 

820 

1400 

1620 

Support 

210 

210 

230 

230 

260 

Fairing 

350 

360 

370 

380 

500 

Actuation 

450 

450 

450 

450 

450 

Total 

1610 

1770 

1870 

2460 

2830 


Table 9.1 Part Count of Trailing-Edge Flaps with Hooked Track Supports 


Part Count 
Component 

Single-Slotted 

Fixed 

Vane/Main 

Articulating 

Vane/Main 

Main/Aft 

Triple-Slotted 

Flap Panel 

600 

750 

820 

1400 

1620 

Support 

150 


164 


185 

Fairing 

100 

100 

105 

110 

145 

Actuation 

300 

300 

300 

300 

300 

Total 

1150 

1300 

1389 

1974 

2250 


Table 9.2 Part Count of Trailing-Edge Flaps with Link/Track Supports 


Part Count 
Component 

Single-Slotted 

Fixed 

Vane/Main 

Articulating 

Vane/Main 

Main/Aft 

Flap Panel 

590 

740 

810 

1380 

Support 

200 

200 

220 

220 

Fairing 

200 

205 

210 

220 

Actuation 

200 

200 

200 

200 

Total 

1190 

1345 

1440 

2020 


Table 9.3 Part Count of Trailing-Edge Flaps with External Hinge Supports 
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Part Count 
Component 

Slat with 
Slave Tracks 

Slat without 
Slave Tracks 

Fixed Camber 
Kruger 

Variable Camber 
Kruger 

Movable L. E. 

1000 

840 

800 

1500 

Actuation 

300 

300 

300 

400 

Fixed L. E. 

1400 

1100 

1000 

1200 

Total 

2700 

2240 

2100 

3100 


Table 9.4 Part Count of Leading-Edge Devices 
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10.0 Sample Application 

The design methodology previously described was applied to evaluate a high-lift 
configuration on a DC-9 type wing. The basic geometry of this test case consists of a 
1 ,000 ft 2 wing with an aspect ratio of 8.5, taper ratio of 0.25, wash-out of 3°, and a 
quarter-chord sweep angle of 24°. The wing thickness varies from 13% at the root to 
1 0% at the tip. The high-lift system consists of a full span slat with a chord of 0. 1 5c and 
a single-slotted flap with hooked track supports, a chord of 0.36c, and a semi-span of 
0.67 b/2 . The aerodynamic performance of the wing with extended trailing-edge flaps is 
evaluated with the slat both retracted and deployed, and system weight and cost are 
determined for both leading and trailing-devices. 

Given: 


S = 1,000 ft 2 

> 

& 

II 

OO 

Ln 

X = 0.25 

s w = 3.0° 


s i 

A = 24° 

— = 0.08 

* 

c 

— =0.36 

Cl c 

— = 0.145 

c 

c 

^TEinb'd = °- 10 

^TEoutb'd =0/77 

^LEjnb'd = 010 

Worth'd =0 - 97 

S TC = 30° = 0.524 rad 

5 le = 30° = 0.524 rad 



Assume for simplicity that the aerodynamic characteristics of the cruise airfoil sections 
remain constant along the span: 
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C, = 5.49 rad-' 

a c!ean 

C| 

°clean 

0.225 

C d . = 0.0067 

mm clean 

c, . 

■min c | ean 

= 0.225 

k 0 = 0.0038 

c, 

' max clean 

= 1.575 

C m =-0.06 

ra °clean 

( d Cml 

— — =-0.039 

^ dC ' ' clean 


Airfoil with extended trailing-edge flat: 


Calculate lift curve slope using the default Fowler motion function: 


f = sin[2( 0.524)] = 0.866 

from eqtn. (18) 

— = 1 + 0.08(0.866) = 1.07 
c 

from eqtn. (17) 

f e . = 1-0.55(0.36) = 0.802 

L I 

from eqtn. (12) 

t - 0.802 

X = j 1 + 0.2 tan"' [- 5.2(0.524) 5 JJ = 0.967 

from eqtn. (11) 

C, a = 0.967(5.49)(1.07) = 5.680 rad' 1 

from eqtn. (23) 


Next, the angle of attack at zero lift is calculated at the flap panel edges for example, but 
a 0 should be calculated at each control point along the wing in order to determine the 
spanwise lift distribution: 
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5 theory 


•-) 

= I V0.36(l - 0.36) + sin _l V036 =0.715 

71 *■ 


fvisc = 0.785 


at inboard section of flap panel, r| TEjnb , d - 0.10 : 


t/c = 0.13 

f t/c = 1 + 1.04(0.13) = 1.135 

a 8 = -0.36 + 0.967[(0.71 5)(1.13 5)(0.785) + 0.36] = 0.604 
AC, te = 5.49(0.604)(0.524)(1.07) = 1.859 
C lo =0.225(1.07) + 1.859 = 2.100 
2.100 

a « = -Iii5 = ' 0 ' 370 rad 

at outboard section of flap panel, t| TEou , b , d = 0.77 : 


t/c = t/c 


root 


^TEputb'd ^ root 
^ — 'H root 




= 0.13 


0.77-0.10 


1 - 0.10 


(0.13-0.10) = 0.108 


f t/c = 1 + 1.04(0.108) = 1.112 

a 5 = -0.36 + 0.967[(0.71 5)(1.1 12)(0.785) + 0.36] = 0.592 
AC 1te = 5.49(0.592)(0.524)(1.07) = 1.822 
C, o = 0.225(1.07) + 1.822 = 2.063 


a 


o 


2.063 

5.680 


= -0.363 rad 


from eqtn. ( 1 0) 
from eqtn. (14) 

from eqtn. (13) 
from eqtn. (9) 
from eqtn. (8) 
from eqtn. (25) 

from eqtn. (24) 


from eqtn. (13) 
from eqtn. (9) 
from eqtn. (8) 
from eqtn. (25) 


from eqtn. (24) 
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Calculate drag polar: 

AC <. mi n TE = 0.038(0.524)' 74 (^)'\l.07) 


= 0.0170 


AC, . 


min te 


2.3(0.524)- 2(0.524) 4 


0.36 

0.3 


- 0.17 


(1.07) 


= 1.09 


Ak p = 


0.00012(0.524)+ 0.0097(0.524) 2 


(0.36) 255 ( 1 


V 0.3 ) V 1.07 


= 0.00406 


C dmjn = 0.0067(1.07)+ 0.01 70 = 0.0242 
C, m . n =0.225(1.07)+ 1.09 = 1.33 

k = 0.00381 — I +0.00406 = 0.00761 
p V1.07/ 

C d = 0.0242 + 0.0076 l(C, - 1.3 3) 2 


Calculate pitching moment about the aerodynamic center: 


x 

— = 0.5- 0.1 34(0.524) 046 = 0.400 
c' 

at inboard section of flap panel: 

AC m = 1.859[0.25 - 0.400(1.07)] = -0.33 1 

TE 

C m = -0.06(1.07) 2 - 0.33 1 = -0.400 

IIIq 


from eqtn. (26) 
from eqtn. (27) 

from eqtn. (28) 

from eqtn. (35) 
from eqtn. (36) 

from eqtn. (37) 
from eqtn. (38) 

from eqtn. (41) 

from eqtn. (40) 
from eqtn. (39) 


at outboard section of flap panel: 
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AC m = 1.822[0.25 - 0.400(1.07)] = -0.324 from eqtn. (40) 

TE 

C mo = -0.06(1.07) 2 - 0.324 = -0.393 from eqtn. (39) 


C, and a 0 are given as input to the Weissinger method in order to calculate the lift 
distribution of the wing. The lift distribution can now be used to calculate C L , C D ., 

C D , and C M . The total aerodynamic coefficients of the wing are given in table 10.1. 

Calculate maximum sectional lift: 

/ rv o 0.16 

f Cf = = 1.030 from eqtn. (63) 

C '8max = 8 ‘ 8 " 7 - 8 (°- 524 ) 0 65 ]( L03 °) = 1786 fr0m ec l tn - ( 62 ) 

AC - = 3.786(0.524X1.07) = 2.1 23 from eqtn. (6 1 ) 

max 

C, max = 1.575(1.07)+ 2.123 = 3.81 from eqtn. (58) 



Calculate lift curve slope using the default Fowler motion function: 


f = sin[2(0.524)] = 0.866 

from eqtn. (18) 

d = sin[2(0.524) 09 ] = 0.899 

from eqtn. (19) 


— = 1 + 0.866(0.08) + 0.899(0.145) = 1.20 
c 


from eqtn. ( 1 7) 
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f and x do not change with the addition of a leading-edge device. 


C, = 0.967(5.49)(1.20) = 6.371 rad' 


from eqtn. (23) 


Calculate the angle of attack at zero lift at the flap panel edges: 

The lift-effectiveness of the flap remains the same with the addition of a slat. 


at inboard section of flap panel, r| TE . = 0.10 : 


AC 1tc = 5.49(0.604)(0.524)(1.20) = 2.085 


from eqtn. (8) 


C, =0.225(1.20) + 2.085 = 2.355 


from eqtn. (25) 


T355 rad 

0 6.371 


from eqtn. (24) 


at outboard section of flap panel, rj TE .,. = 0.77 : 


AC Itc = 5.49(0.592)(0.524)(1.20) = 2.044 


from eqtn. (8) 


C, =0.225(1.20) + 2.044 = 2.3 14 


from eqtn. (25) 


_ 2314_ rad 
0 6.371 


from eqtn. (24) 


Calculate drag polar: 


0 36^ 1-4 

AC h . =0.038(0.524)' 74 — — (1.20) = 0.0191 from eqtn. (26) 

a min-rr: x \ f\ ) 


AC - minTE =[2.3(0.524)-2(0.524) 4 ][— J (120) from eqtn . (27) 


= 1.23 


2 Y 0.36^ 2 ‘ 55 7 1 


Ak p = 0.00012(0524) + 0.0097(0.524) I— I I ttt! from eqtn. (28) 


= 0.00362 



76 


AC, . = 0.0074(1.20) = 0.0089 

a min LE 

from eqtn. (32) 

AC, . =0.76(1.20) =0.91 

'mmyj? 

from eqtn. (33) 

/ 

Ak n =0.00381 

P LE V 

-M =0.00318 
, 1.20/ 

from eqtn. (34) 

C, . =0.0067(1.20) + 0.0 191 + 0.0089 = 0.0360 

a mm 

from eqtn. (35) 

C, . =0.225(1.20) + 1.23 + 0.91 = 2.41 

'min 

from eqtn. (36) 

k P=°- 0038 (,.20. 

) + 0.00362 + 0.003 1 8 = 0.00997 

from eqtn. (37) 

= 0.0360 + 0.00997(C, -2.4lf 

from eqtn. (38) 


Calculate pitching moment about the aerodynamic center: 


— = 0.400 
c' 

from eqtn. (41) 

0 LE = cos -1 [1-2(0.145)] = 0.781 rad 

from eqtn. (44) 

C m - — sin(0.78l)[l -cos(0.78l)] = -0.102 rad' 1 

5lE 2 

from eqtn. (43) 

AC m = -0.1 02(0.524)( 120) 2 = -0.0770 

0 r c 

from eqtn. (42) 


at inboard section of flap panel: 

AC m = 2.085(0.25 - 0.400(120)] = -0.480 from eqtn. (40) 

TE 


C m = -0.06(1.20) 2 -0.480-0.0770 = -0.643 
IlI 0 


from eqtn. (39) 
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at outboard section of flap panel: 


AC m = 2.044[0.25 - 0.400(1.20)] = -0.470 

TE 

from eqtn. (40) 

C m = -0.06(1.20) 2 - 0.470 - 0.0770 = -0.633 

1 "o 

from eqtn. (39) 

Calculate maximum sectional lift: 


AC, =3.786(0.524X1.20)= 2.381 

max yg v 

from eqtn. (61) 

Aa max = 0.023 + 1.7(0.145)- 1.8(0.145) 2 -7.3(0.145) 3 
= 0.209 

from eqtn. (60) 

AC, =6.371(0.209)= 1.332 

'max LE 

from eqtn. (59) 

C, =1.575(1.20)+ 2.381 + 1.332 = 5.60 

‘max v 7 

from eqtn. (58) 

This information is insufficient to determine the optimum high-lift system, so the 


weight and cost data is calculated next. 


Calculate weight and cost of trailing-edge flap: 


C TE inb . d ~ 0.36c jnb'd 
= 0.36c 


= 0.36- 


root 

2 ' 


1 + X V AR 


= 0.36 


1000 


1 + 0.25 V 8.5 


= 6.248 ft 
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C TE 0Utb 'd -° 36c outb'd 

= 0.36c 
= 0.36 


1 ^ TEputb'd ft root ^ ^ ^ ^ 


1-T) 


root 


= 0.36- 


1+X VAR 
2 


j ^TEputb'd ^Iroot ^ ^ 


1000 


1 + 0.25 V 8.5 
= 2.759 ft 


1 — 


1 ft root 

0.77-0.1 

1-0.1 


(1-0.25) 


S te -(c TEinb , d +C TE 0Ut b'd j^TE 2 

_/ \ VSAR 

\ TEjnb'd ^TE 0U (b'd jftTE 2 

, . J 1000(8.5) 

= (6.248 + 2.759)0.67 

= 278.2 ft 2 

Since f fow =l, equation (72) reduces to: 

W TE =(a, +a 2 +a 3 + a 4 )S XE 
W TE = 8.9(278.2) = 2476 lb 
PC XE =1610 

% 

Cost XE = 1.888 1(2476)(1 6 1 0) 07 = $ 0.821 x 10 6 
Calculate weight and cost of slat: 


c LE inb 'd - 0.145cj n b'd — 0.145c root 


= 0.145 



= 0.145 


1 + 0.25 




iuUU 

8.5 


= 2.516 ft 


from eqtn. (80) 
from eqtn. (79) 
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c LE 0U tb'd 0-145 c oll tb'd 
= 0.145c 


' root 


^ ^ LEputb d ^root ^ ^ 


i-n 


root 


= 0.145- 


= 0.145- 


\+X V AR 
2 


^ ^ LEpmb’d ^ root ^ ^ 


i-o 


root 


liooo 

r 0.97-0.1. ,1 

1- , —(1-0.25 

8.5 

L i-0-i J 


= 0.692 ft 


S LE “ ( C LE inb'd + C LE 0 utb'd LE 2 


\ VSAR 

^EEjnb'd ^ LE 0U tb'd JOlE 2 


, , Jl 000(8.5) 

= (2.5 1 6 + 0.692)0.87 j- 

= 128.7 ft 2 

Equation (78) can be expressed as: 


W LE — ( a 5 +a 6 +a ?)SLE 


W LE = 15.3(128.7) = 1969 lb 


PC LE =3100 

from eqtn. (81) 

Cost LE =1.7339(1969X3100)° 7 =$ 0.949 xlO 6 

from eqtn. (79) 


A summary of the results is shown in tables 10.1 and 10.2 and figures 10.1 and 10.2 
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Single-Slotted 

Single-Slotted with 
3 Position Slat 

c c 

1.28 

1.38 

Co (a=0°) 

0.0953 

0.1228 

C L 

max 

1.95 

3.03 

^max’ 

9.32 

21.1 

Table 10.1 Aerodynamic Performance ol 

'a Single-Slotted Flap 


With a Track Support (With and Without a Slat) 



Single-Slotted 

3 Position Slat 

Total 

System Weight, lb 

2476 

1969 

4445 

System Cost, $ 

0.821 x 10 6 

0.949 x 10 6 

1.770 x 10 6 


Table 10.2 System Weight and Cost Data of a Single-Slotted Flap With a Track Support 

and a Slat With Slave Tracks 
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Figure 10.2 Drag Polars of Sample Test Cases 
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11.0 Concluding Remarks 

Today’s highly competitive market necessitates the ability to develop a high-lift 
system that allows a civil transport aircraft to meet the stringent FAR Part 25 takeoff and 
landing requirements with minimal direct operational cost. This presents a daunting task 
for designers, since the optimization of a high-lift system is extremely complex. This 
report presents a comprehensive design methodology which includes aerodynamic 
performance, structural weight considerations, system complexity, and cost to be used 
within an aircraft conceptual design tool such as ACSYNT. Thus, aircraft performance 
and cost can be analyzed in parallel to select the optimal configuration. 

Similar methodologies have been created in the past, but they should not be used 
today for two reasons. The methods generally consist of outdated experimental data, and 
the data is based on low Reynolds numbers. Today, CFD offers the best alternative, 
since it would be too costly to generate a large high-Reynolds-number experimental 
database as was done in the past. Computational methods can accurately, yet 
inexpensively, predict the flow field around multi-element airfoils in a reasonably short 
period of time. However, the most important advantage of using CFD is that airfoils can 
be analyzed at realistic flight Reynolds numbers. 

CFD was the primary tool used in this research project to construct a database 
consisting of the NLR-7301 two-element and the Douglas LB-546 three-element airfoils. 
However, supplemental configurations such as a double-slotted and triple-slotted flap 
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must be added to build a more comprehensive database. This can be done without 
difficulty, since the new design methodology is easily extensible. 

From this database empirical equations were constructed which are used to 
calculate two-dimensional aerodynamic coefficients. A spanwise lift distribution is then 
determined from a modified lifting-line theory using the two-dimensional data as input. 
Thus the total lift, drag, pitching moment, and maximum lift of the wing can be computed 
faster than with a panel method while incorporating the viscous effects of slotted high-lift 
devices. 

Since the design methodology was intended to be used in an optimization routine, 
a reasonable compromise between CPU requirements and the accuracy of the predictions 
had to be found. Consequently, some assumptions and simplifications were made. One 
such simplification rests in the prediction of weight using historical data. Although this 
technique offers a consistent approach for predicting weight, it is a very simple method. 
Further research should be conducted to consider a method which incorporates loading 
information from CFD with slender beam analysis to determine the optimal weight of a 
configuration. 

This methodology also neglects spanwise viscous flow effects since it employs a 
modified lifting-line method. Clearly, this problem cannot be solved until full 3-D 
viscous flow solutions become practical. Using the three-dimensional Reynolds-averaged 
Navier-Stokes (RaNS) equations will not only account for spanwise effects, but will also 
allow for inverse wing design. Nevertheless, until the 3-D RaNS equations for multi- 
element wings at high-lift conditions can be solved in a reasonable time period (i.e. less 
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than a day) the methodology presented in this paper will be a useful tool in the 
development of high-lift systems on civil transport aircraft. 
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Appendix A 

Spanwise Load Distribution Calculations of a Wing with 
Arbitrary Plan Form Using the Weissinger Method 


Summary 

The application of Weissinger's theory is shown for calculating the spanwise 
loading across a wing. From this procedure a spanwise load distribution can be 
calculated for a wing of arbitrary plan form from a 2-D lift curve slope and angle of zero 
lift. The geometry of the wing may be asymmetric about the centerline and can include 
sweep, taper, or twist. Partial span flaps may be modeled as discontinuities in twist. 
Fowler motion can also be modeled as an increase in chord length. Lift, induced drag, 
pitching moment, and induced angle of attack may then be calculated. 

Method 

The Weissinger method models the wing as a plate of zero thickness, but the 
planform and twist remain identical to the actual wing. The chordwise load distribution 
at each span station is concentrated into a lifting line located at the wing's quarter-chord 
line. The theory requires that the quarter chord line is straight, however a discontinuity is 
allowed at the plane of symmetry so that a swept wing may be considered. 

The method assumes that the lifting line and its trailing vortex sheet are 
continuous. However, discrete values for the circulation strength are determined only at 
the span stations along the quarter-chord line corresponding with control points. The 
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number of control points, m, is specified, and each is placed at the three-quarter-chord 
line. Mathematically this implies that the lift curve slope is 2n rad' 1 , and a correction 
method must be incorporated for a lift curve slope that varies from the theoretical value. 

The boundary conditions specify that the induced angle of attack due to the 
downwash of the trailing vortex is equal to the angle of attack of the plate. This ensures 
tangency of the flow to the plate at each control point. 

Substituting the wing with a bound vortex system and applying the boundary 
conditions enables the formation of a set of simultaneous equations. Each of these 

C,c 

equations is a function of the angle of attack of the wing, the load coefficient G ■ = — — 

3 2b 

at each span station j on the quarter-chord line, and the influence coefficients Ay which 
are only a function of geometry and relate the influence of the circulation at any point j 
along the lifting line to the downwash at any control point i. The distribution of the load 
Gj may then be calculated from the set of m simultaneous equations. 

m 

a i=X A i,j G j’ i = 1,2,3, ...m (A-l) 

j=i 

where 


Wj 

— -y — (^wing 

’ on 


(A-2) 


Each equation gives the downwash angle at the control point r| = cos<j> j resulting from 
the circulation effects of m points along the span. 

The following trigonometric substitutions are used to define span locations at the 
control points or on the lifting line respectively. 
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(A-3) 

(A-4) 

iz:cients are shown below. 


(A-5) 

(A-6) 



(A-7) 


- 3 i 

— -os^) 2 2(m+l) 


ents are given as 


f + 1 f M 

-i.oM.o + H.M+rj.M+l v -1 , . 

o + 2. L i,kfj,k) 

z k=l 


(A-8) 


(A-9) 


— of span stations used to integrate the down wash and does not have 

zz z i The subscripts 0 and M+l denote values for k. Also note that 

_ = j (A- 10) 


in coefficients L ik are defined as follows: 


( a in 

V A A 
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k7t 


n = 


COS 

m+ 1 


(A- 12) 


for r\ < 0 


L..k = 


1 V[l + (b/c)(|r|l + r|)tanA] 2 + (b / c) 2 (r| - r|) 2 

(b/c)(r|-r|) l+(b/c)(|r|| + r|)tanA 


+ 


2tanA-J[l + (b/c)|ri|tanA] 2 +(b/c) 2 r| 


(A-13) 


[1 + (b/c)(|ri|- r()tanA][l + (b/c)(|ri| + ri)tanA] 


for r| > 0 


kk = 


1 fN /[l + (b/c)(|r||-ri)tanA] 2 + (b / c) 2 (r\ - r|) 2 


(b / c)(r| — n) ^ 


1 + (b / c)(|ri| - r|)tan A 


1 } 


(A- 14) 

The mathematical series coefficients f j k are used in a numerical integration method and 
are shown below: 


(|) k = cos 


krc 

M + l 


(A-15) 


2 m 

fj,k =^— jX^sin^jCos^k 


(A- 16) 


Sectional Lift curve Slope Correction 

The Weissinger method assumes that the lift curve slope is 2n rad - 1 . This implies 
that the control points lie along the three-quarter-chord line. When considering airfoils 
with lift curve slopes that vary from this theoretical value (as a result of viscous, 
compressibility, or geometry effects), a correction method must be incorporated. This 
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involves moving the control points from the three-quarter-chord line; but before this is 
done a derivation of this location will be shown. 

The velocity induced by an infinite vortex of strength T and distance h is given as 


w = 


2 7th 


(A- 17) 


and the sectional lift coefficient is related to the circulation strength as follows: 


pV oo r=^pVicC, (A- 18) 

Rearranging these two equations and combining them yields 
w h 

C,=47t- — (A-19) 

V. c 


w 


Since a = the above equation can be rewritten as 


C, h 

— = 4tt- 
a c 


(A-20) 


h 1 

If the lift curve slope is assumed to be 27t then it can be shown that — must be — c . Then 
F c 2 


the point at which no flow passes through the plate occurs at -c from the lifting line or 

•» 

at the three-quarter-chord line. 

If the lift curve slope deviates from the theoretical value each control point 
location must be modified. This means that each control point must be moved forward of 
the three-quarter-chord line if the slope is less than 2n rad" 1 or must be moved rearward if 

b 

it is greater. This variation can also be modeled as a change in the value of — by the 

c v 
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a eX p 

ratio of the experimental lift curve slope to the theoretical lift curve slope — — at the 

desired Mach number, M, where P = y[\ - M 2 . The modified geometry can be expressed 
as 



(A-21) 


where 


k v 


2n/p 


(A-22) 


Caution should be taken when large deviations in the lift curve slope are used. The 
method does not rigorously allow such modifications, and the accuracy of the results may 
deteriorate for large angles of sweep. 


Compressibility Correction 

Compressibility effects are included with the aid of the Prandtl-Glauert rule. It 
approximates the effects of compressibility as an increase in the local chord and sweep 


1 

angle by the factor — . 



J__b 

P c v 


(A-23) 


Ap = tan 1 ( 


tan A 



(A-24) 


Thus the span loading for the original plan form at a given Mach number can be 


calculated for an incompressible flow using a modified geometry. 
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Aerodynamic Characteristics 

The circulation strength is calculated directly from the dimensionless load 
coefficient and has the dimensions [ft 2 /s]. 

r v =bV 00 G v (A-25) 

The sectional lift coefficient is also calculated from the load coefficient as 


b 

C , v = 2 -Gv 

'■'V 

The total lift coefficient is determined from 
ttAR ^ 

C L =— rZ. G v sin <() v 

m + 1 V= | 

The induced angle of attack at the quarter-chord point of station v is given by 


(A-26) 


(A-27) 


^iv b vv G v b vn G n 


(A-28) 


n=l 


where the primed summation sign indicates that the value for n = v is not summed. The 
total induced drag is calculated from 
7tAR ^ 

C Dj = — 7 TZ. G v«i v sin(t) v 

m T I v _| 


(A-29) 
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Appendix B 

Weight Estimation of Other Control Surfaces 

Like leading and trailing-edge device weights, the estimation of aileron and 
spoiler weights is based on historical data of existing aircraft. Similarly, the total system 
weights are broken up into components of surface, support, and actuation. The weight 
estimation equations given below are only a function of the planform area of the 
representative control surface. 

Aileron Weight 


Surface Weight: 

w s „„ =5.5S,„ (B-l) 

Support Weight: 

W„ w =0.74S u , (B-2) 

* 

Actuation Weight: 

W„=4.3.S,„ (B-3) 

Total Aileron Weight: 

W,i, = w surf + W sppp + W» (B-4) 
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Spoiler Weight 


Surface Weight: 

W 5u rf =110 + 1.65 - S sp (B-5) 

Support Weight: 

W^ up p = 1-3 ■ S sp (B-6) 

Actuation Weight: 

W act = 270+ 1.35 -S sp (B-7) 

Total Spoiler Weight: 

W sp =W sm(+ W SUH , + W„ (B-8) 
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Appendix C 

High-Lift Weight Data of Various Aircraft 

(All weights are in pounds.) 


Trailing-Edge Device Weight Data 





Spoiler Weight Data 
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Appendix D 


Database Formulation 


The first step in integrating a database into the high-lift module is to formulate the 
data based on important parameters. This has already been done, and the empirical 
equations can be found in the documentation. The second step is to determine 
appropriate empirical coefficients for the formulas. This has been completed for leading- 
edge devices and single-slotted trailing-edge flaps but can be done for other 
configurations by finding the empirical coefficients which minimize the error between the 
predictions and the actual data. In other words, iterate on the empirical coefficients to 
minimize the RMS: 

RMS = yaverage(predicted values - database values) 2 

The empirical equations for x, AC dmin , AC W Ak„, C, s , and x cp are given in equations 1 1 , 

26, 27, 28, 62, and 41 respectively and must be fit to a new configuration data set in order 

to integrate the new configuration. The database values of these parameters are 

% 

determined from the CFD data as follows: 


a 


8 


AC 


>TE 


C 

5 xe — 

'“clean 1 fc c 


a 8 -) 




(Find x which accurately predicts 



